Single-atom laser generates nonlinear coherent states 
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The stationary state of a single-atom (single-qubit) laser is shown to be a phase-averaged nonlinear 
coherent state — an eigenstate of specific deformed annihilation operator. The solution found for 
the stationary state is unique and valid for all regimes of the single-qubit laser operation. We have 
found the parametrization of the deformed annihilation operator which provides superconvergence 
in finding the stationary state by iteration. It is also shown that, contrary to the case of the usual 
laser with constant Einstein coefficients describing transition probabilities, for the single- atom laser 
the interaction-induced transition probabilities effectively depend on the field intensity. 

PACS numbers: 42.50.-p, 32.80.-t, 42.55.-f, 42.60.Da 
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I. INTRODUCTION 

Recently the 50th anniversary of the invention of the 
laser [l[ has been celebrated, while the ideas lying behind 
light amplification [2j are almost 100 years old. In the 
beginning of the laser era it was realized that laser pho- 
tons are emitted in specific coherent superpositions - 
coherent states |3J , which form a kind of border between 
classical and nonclassical states of light. 

The general tendency of miniaturization of electronic 
and optical devices and components is also observable 
in the diminishing of the laser size down to the value 
of the wavelength. The use of microcavities of differ- 
ent types, like interferometric and Fabry-Perot microcav- 
ities, microcolumns, whispering gallery mode resonators 
(microdisks and microspheres), 2D- and ID-tapered pho- 
tonic crystal resonators, allows the single-mode thresh- 
oldless regime of lasing to be reached due to the increase 
of the ratio of photons spontaneously emitted into the 
lasing mode to the number of photons emitted into non- 
lasing modes. The extreme case of the active element of a 
microlaser (or micromaser) is a single emitter — an atom 
(in. Rydberg 0, Q or lower electronic states [H, 0] ) , an ion 
Q, a quantum dot 0] or a superconducting qubit, play- 
ing the role of an artificial atom in an electrical resonator 
in recent demonstrations of a single-qubit laser [Icj |. 

The one-atom-one-mode microlaser is of great impor- 
tance as a limiting case of lasers. This intrinsically 
quantum system with a number of properties very dif- 
ferent from those of ordinary lasers requires specific cav- 
ity quantum electrodynamics methods for its description 
[ill - tl3l | . Rabi splitting [l4], HH , the collapse-and- revival 
phenomenon [lg, [TtJ , and the photon blockade effect [l8[ 
are a few examples of quantum effects observed in the 
system (see also fl9l|). 

In contrast to conventional lasers, microlasers (and es- 
pecially single-atom lasers) are known to be sources of 
nonclassical light [20I - I23I ] . It has already been shown that 
a single- atom laser, considered within the scope of the 
strong-coupling regime, can produce special kind of non- 
linear coherent states (NCSs), namely, Mittag-Leffler co- 
herent states [24| • In this paper we provide a general uni- 
formly applicable description of the single-atom laser and 



show that it generates NCSs for any values of the interac- 
tion parameters. A NCS can be written as an eigenstate 
of a specific deformed annihilation operator. It should be 
emphasized that the solution found is unique and follows 
from the master equation exactly, without any approxi- 
mations. We believe that the finding is both interesting 
from the fundamental point of view (as a connection be- 
tween the classes of deformed annihilation operators and 
a single-qubit laser) , and useful for further analytical and 
numerical investigations of stationary state non-classical 
properties, not accounted for correctly by approximate 
solutions. 

In the case of strong coupling our solution agrees with 
the corresponding approximate solutions [H, [25[, pre- 
dicting, however, state nonclassicality, not described cor- 
rectly by the strong-coupling approximation, in regimes 
of weaker coupling. It is worth noting that, although 
nonlinear properties of (multi-emitter) lasers have been 
investigated for quite a long time (see, e.g., Refs. [26l - [28| ). 
the nonlinearity and nonclassicality of the properties con- 
sidered here are new and characteristic of the inherently 
quantum nature of a single-emitter laser. 

The intrinsic quantum character of the light-matter in- 
teraction in single-atom lasers reveals itself in the impos- 
sibility of describing the lasing effect be means of field- 
independent spontaneous and induced transition proba- 
bilities, as in the case of a conventional laser. The ef- 
fect has been mentioned for the strong-coupling regime 
in Ref. [Hj]. Here we show that this property is general 
and is preserved also beyond the strong-coupling regime. 
We present both numerical and uniformly applicable an- 
alytical expressions for the transition probabilities which 
are intensity dependent and provide an explanation of 
the found "saturation" effect. The observed features of 
a single-atom laser are a manifestation of its quantum- 
ness, revealing itself in an extremely strong correlation of 
atom and field states (compared to conventional lasers) 
and leading to invalidity of mean-field and other semi- 
classical approaches. 

The paper is organized as follows. First, we intro- 
duce the model of an incoherently-pumped single-atom 
laser and derive the equations describing its stationary 
state. Then an analytical solution of the equations in 
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the form of generalized coherent states is provided. The 
solution is obtained by introducing the state-dependent 
operator d(n), which describes the difference between 
the exact solution and the solution obtained under the 
strong-coupling approximation. It is shown that the it- 
eration scheme for finding d(n) is unconditionally stable. 
Moreover, the scheme does not depend on the bound- 
ary values of d(n). We demonstrate the use of the it- 
eration method both for numerical calculations and for 
constructing uniformly applicable analytical approxima- 
tions. Then, specific properties of the stationary-state 
nonclassicality are discussed on the basis of phase space 
quasi-distributions. In the last section we discuss the in- 
terpretation of the system evolution equations in terms 
of spontaneous and induced transition probabilities. We 
show that for a single-atom-single-mode system the in- 
tracavity spontaneous emission probabilities strongly de- 
pend on the number of photons in the mode (in contrast 
to usual case, when the normalized probabilities are con- 
stant), which is a manifestation of the inherently quan- 
tum features of single objects. 

II. EQUATIONS 

A single-atom laser is considered within the framework 
of a model system consisting of a two-level atom with the 
ground state |1) and excited state |2), interacting with 
a resonance field mode with coupling constant g. The 
atom is pumped incoherently with mean rate R\ 2 - In 
addition, decay of the resonance field mode and decay 
and dephasing of the atom with rates k, R21, and T, 
respectively, are taken into account. 

The master equation for the density matrix, reduced 
over the states of the environment, in the interaction rep- 
resentation has the form: 

i 

P= -t [H,p\ + 2nL a p + Ri 2 L a p + R 21 L a ._p + TL a . z p, 
n 

(!) 

where the operators a +1 <r_, a z and <v , a describe the 
dynamics of the atom and the field, respectively, and the 
relaxation is described by Lindblad operators: 2Lxp = 
2XpX^ - X^Xp - pX^X. The atom-field interaction is 
described by the Jaynes-Cummings Hamiltonian: H = 
Kg (a T er_ + acr + ) . 

In the paper we investigate the properties of the sta- 
tionary state of the system. The following four nor- 
malized parameters are used below for simplifying the 
equations: — R\2/(Ak), Vq — {R21 — 2k)/ (4k), /io = 
o,q + uq + T/k and r\ = g 2 /n 2 , describing the pump, the 
atomic loss excess over the field loss, the dephasing and 
the atom-field coupling, respectively. 

Introducing the jump (J) and photon number (N) su- 
peroperators 

J p = apa\ N p = aJap 
and decomposing the density matrix in terms of atom 



states as 

p = pn®|l )(1|+/J22®|2 )(2\+a(v + iu)®\2)(l\+h.c, 

one can find the following properties of the stationary 
state. 

a. The stationary state is unique. This statement 
follows directly from the form of master equation and 
can be proved by considering the evolution of the trace 
distance, defined as D(p, a) — i Tr (\p — a\), where \A\ — 

v 7 A^A (see, e.g., Ref. [Hj]). It is known that the station- 
ary state is unique, when the condition D(p(t),a(t)) < 
holds for any non-equal solutions p(t) and a(t) of the 
master equation. 

Eq. consists of two parts: the Hamiltonian one and 
the sum of Lindblad form superoperators. The Hamilto- 
nian part does not influence the distance between quan- 
tum states and, therefore, preserves the uniqueness prop- 
erty of the stationary state. The Lindblad part of the 
master equation does not include atom-field interaction 
and describes independent interactions of the atom and 
the field with thermal baths. The evolution, cause by 
these interactions only, has the unique stationary state 

po = |0)(0|® {R 2 i\l){l\+ Ri 2 \2){2\) I {R 12 + R 21 ) 

and corresponds to a strictly negative time derivative of 
the trace norm: 

D(p(t),a(t))<Ohxp(t)^a(t). 

Finally, the above condition is satisfied for the evolution, 
described by Eq. ^ , and the stationary state of a single- 
atom laser is unique. 

b. The operators pn, p 22 , u and v, acting on the field 
mode, are diagonal in Fock basis. Diagonality of the 
operators follows from the stationary state uniqueness. 
Eq. ([T]) is invariant under the transformation a —> ae~ 4< ^, 
at ->■ ate**, er± -> a ± -e ± ^. Therefore, if the start- 
ing state is described by the diagonal operators pn, p 22 , 
u and v, the stationary state will also possess the diag- 
onality property. Together with the uniqueness of the 
stationary state, it implies that the operators pn, p 22l 
u and v become diagonal in the limit t — ► 00 for any 
starting state. 

c. The operator v vanishes in the stationary state. 
This operator evolves independently of operators the pn, 
p 22l and u: 

v = 2nL a v - |k + i (i?i2 + #21) j v. 

The first term of the equation describes trace-preserving 
dissipative dynamics, the second one corresponds to ex- 
ponential decay of the operator v. 

d. The operator u is defined in the a unique way by 
the field density operator Pf = pu + p 22 : 

u = Pf/Vv- 
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FIG. 1: Scheme of energy levels and transitions: wavy arrows 

— mode decay, double arrows — atomic excited state decay, 
single arrows — pumping. The balance of the transitions, 
crossing dashed lines, is described by Eqs. and (|20p : (a) 

— transitions between states with n and n + 1 excitations 
(atom + field); (b) — transitions between states with n and 
n + 1 photons. 



This property follows from the master equation and the 
operators diagonality in the stationary state. 

The operators pu and P22 satisfy the following equa- 
tions: 

{2v + N + l) p 22 = (2flo - J) Pu, (2) 



(N + l)p 2 2 = Jjpii + -(Mo + -W- J)p f >. (3) 

If Eq. © has a simple interpretation as the balance of 
the number of total excitations in the system [Fig. QJa)], 
Eq. ([3]) has a more complex interpretation and can be 
considered in terms of held-induced transitions between 
ground and excited states of the atom. In the limiting 
case of weak atom-field correlation the transitions corre- 
spond to ordinary spontaneous and induced transitions. 

Because of the strict positiveness of all elements 
Pu{n) = {n \pii\n), following from Eq. (JlJ, it is possi- 
ble to define a superoperator d(N ) , diagonal in the Fock- 
state basis, by the following equation: 

d(N)p n = -(jM> + N-J)p f . (4) 
V 

[It should be noted that for any function f(n), defined 
for n = 0, 1, . . ., the action of the superoperator f(N) on 
diagonal density matrices is also correctly defined. The 
superoperator l/(N + 1) is an example.] 

Using Eq. one can rewrite Eq. §3§ in a simpler way, 

P22 = 1 ^J{l + d(N)}p 11 , (5) 

showing directly that the excited-state photon statis- 
tics p22(n) and shifted ground-state statistics p\\(n + 
1), equalized by frequent intracavity transitions in the 



strong-coupling regime, become unequal in the general 
case: P22(n) / pn(n + 1) = 1 + d(n + 1). It is the function 
d(n + 1) that describes the deviation of the ratio from 
unity. 

III. GENERATION OF NONLINEAR 
COHERENT STATES 

Substituting Eq. <j5j) into Eq. @, we arrive at the fol- 
lowing equation for the conditional density matrix pu- 

A Fll puA\, u = alpu, (6) 

where 




is a deformed annihilation operator [30, 31] with its prop- 
erties completely determined by the discrete function 
F(n) (the deformation function). For the ground-state 
conditional operator pu this function equals 

*li(») = ^+Q + ^){l + d(n)}. (7) 

and is determined by the parameter v>q and the discrete 
function d(n). 

Eigenstates of deformed annihilation operators are 
known as nonlinear coherent states (30L l3lT ] and represent 
a particular case of generalized coherent states (see, e.g., 
|32|). In the special case vq — 0, d(n) = 0, eigenstates of 
the operator A are ordinary coherent states. For vq ^ 0, 
d(n) = (the strong-coup ling regime), the eigenstatates 
are Mittag-Leffler states [24j. In the general case, the 
eigenstate \ao;F), corresponding to an eigenvalue ciq of 
the operator Ap, has the following Fock decomposition: 

oo n n 1 

|oo; F) = const ■ £ |n> ^= J[ -==. (8) 

„=0 m=i V F ( m ) 

It follows from Eqs. ©-© that the density matrix pu 
represents a phase- averaged NCS: 

pu = diag (\a ; F n )(a ;Fn\) . (9) 

Eq. ([5]) implies that conditional (P22) and unconditional 
(pf) field operators also correspond to phase-averaged 
NCSs, but with different deformation functions i 7 22(«) = 
Fu{n)ip{n)/y(n + 1) and F f (n) = F n (n)[l + <p(n)]/[l + 
ip(n + 1)], respectively, where f(n) = Oq / '[Fu(n)n(n)] 
and h(n) = n/[l + d(n)]. 

It is worth noting that the derived representation of the 
stationary state follows from the exact master equation 
([T]) without any additional assumptions and approxima- 
tions and is valid for all values of the system parameters. 
The solution found is general and has the same form for 
all of the five possible regimes of single-qubit laser oper- 
ation: linear, nonlinear quantum, lasing, self-quenching, 
and thermal [33[ (see also Fig. [5] below). 
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IV. CALCULATION ALGORITHM 

Eqs. (H])-© derived above imply that determination of 
the stationary-state density matrix is equivalent to find- 
ing the discrete function d(n). According to Eqs. (J4J)— ([6j) , 
the deviation function d(n) satisfies the following system 
of recurrence equations: 

2 r 

d(n) = - (no + n) {1 + ip(n + 1)} 

11 n (10) 

-n(n + l)ip(n + 1) {1 + tp(n + 2)} 

The value d(n) depends on d(n+l) and d(n + 2), implic- 
itly present in n(n + 1), ip(n + 1) and ip(n + 2). 

Generally, in order to calculate d(n) for n = 0, . . . , no, 
one needs to know correct values of the deviation function 
d(riQ + 1) and d(no + 2) near the starting point n = uq. 
However, the following characteristic properties of the 
map (fl~0)) enable us to calculate the values of d(n) with 
arbitrarily high accuracy without any prior knowledge: 
(i) the map is stable — small deviations from the cor- 
rect solution decrease approximately exponentially dur- 
ing steps, described by Eq. (ITUl) ; (ii) the value d(n), de- 
fined by Eq. (fTUl) . is bounded: 



< d(n) < - (mo + n) 1 + 

r\ \ 2^o + «+l 



2al 



(11) 



This means that the first iteration step brings d[n) close 
to the correct solution regardless of the chosen initial 
values d(n + 1) and d(no + 2), provided these values are 
positive. Fig. [2] illustrates fast convergence of the numer- 
ical solution for d(n) for different initial conditions [the 
bounds shown by the grey region are obtained itcratively 
on the basis of Eq. (fTTj)]. 

The above-discussed stability of the map, defined by 
Eq. (|10[) . provides also a quite simple way to decompose 
d(n) analytically in terms of small parameters by succes- 
sive improvement of the approximations. For example, 
for n>l and all values of the coupling parameter r\ the 
following expression is valid: 



d(n) — — < n 



2n 



n 



oil 

n 



(12) 



For rj >• n Eq. (Ti~2"j) implies that d(n) sa |(n + /Jo), 
which corresponds to the adiabatic approximation [25| . 
This expression, as well as inequality ([TTj) . provides the 
condition for validity of the strong-coupling approxima- 
tion: one can take d(n) ~ for n, *C r\. Fig. [3] shows 
the density matrix elements of the single-qubit laser sta- 
tionary state, calculated on the basis of the above ex- 
pression (dashed, dotted, and dot-dashed lines) and us- 
ing the strong-coupling approximation (solid lines). The 
difference between the exact and approximate solutions 
becomes significant for large photon numbers n and for 
small values of the coupling parameter 77. Fig. 0] shows 
the dependence of the distance between the exact and 
approximate solutions on the coupling parameter r\. 
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FIG. 2: Function d(n): black and grey points — results of 
numerical calculations on the basis of Eq. (|10[) with different 
starting values of n and d(n), d(n+l) (black points — starting 
from n = 40); dashed black line — approximate analytical 
expression [Eq. l|12[l]. valid for n > 5; grey region — values 
of d(n) after the first iterations step [analytical expression, 
improved version of Eq. Parameters: vo = 1, do — 1> 

= 3, 77 = 5. 




FIG. 3: Field density matrix pf. dashed, dot-dashed, and 
dotted lines — results of numerical calculations; points - 
analytical calculation on the basis of Eq. (|12[) for the same 
sets of parameters, as lines; solid lines — calculation on the 
basis of strong-coupling approximation. Parameters: gray 
lines: vq — 0, Oq = 5, fio — 5, r) = 30, 200 (dot-dashed and 
dashed lines, respectively); dark gray lines: vq — 1, ap = 1, 
fio = 3, 77 = 5, 15, 50 (dotted, dot-dashed and dashed lines, 
respectively). 



V. NON-CLASSICAL PROPERTIES OF THE 
STATIONARY STATE 

For the limiting case of highly-excited states (n 3> 77) 
Eq. (JX3J) implies that d(n) ~ 2n/?7 and 



(n I a;Fu), (n \ a;F 22 }, (n | a;F f ) 



a n rj n 



(13) 



5 



1.0 

g 0.6 

3 0.4 
Q 

0.2 




function correspond to s 
tively. 



1, s = and s = — 1, respec- 







20 



40 



60 



80 



FIG. 4: Trace distance between the field density operator 
Pf, calculated numerically, and the operator psc, obtained in 
the strong-coupling approximation. Parameters: gray lines: 
Uo — 0, do = 5; black lines: uo = 1, = 1; solid lines: (iio = 5; 
dashed lines: /Lto = 10. 



In this case the decrease of the density matrix elements 
Pf(n) ~ l/(n!) 2 with growth of n is faster than for any 
ordinary coherent state with nonzero amplitude. This 
fact indicates nonclassicality of the stationary state: any 
classical state can be represented as a mixture of coher- 
ent states with positive weights 0; its matrix elements 
(n \ p\n) cannot decrease faster than for a certain coher- 
ent state with growth of n. It should be noted that the 
strong-coupling approximation predicts decrease of the 
density matrix elements proportionally to psc( n ) ~ 1/nl 
[25}. These are the "tails" of the photon number distri- 
bution, present in the approximate solution and absent 
in the solution found in our paper, that cause a non- 
zero distance between the exact and approximate den- 
sity operators. The numerically calculated trace distance 
D(pf, psc), shown in Fig. |H almost coincides with the 
total weight of the excess " tails" of the approximate so- 
lution. 

To characterize the types of nonclassicality of the sta- 
tionary state, it is useful to consider nonclassicality pa- 
rameters [H], based on considering s-parametrized 
phase-space functions [36T!38j P(a; s), equal to the mean 
value of an observable 



6(a- 



7T(1 



exp 



2(a t -a*)(a-a) 



1 



(14) 

where the colons denote normal ordering of the field op- 
erators. Any of the functions P(a; s) represents a con- 
volution of the Glauber function P(a) with a Gaussian 
weight function: 

(15) 

The Glauber P function itself, Wigner function and Q 



For any classical state the Glauber function is well de- 
fined (except for <5-function type singularities) and takes 
non-negative values. The weight function in Eq. (|15[) 
is strictly positive. Therefore, any classical state is 
characterized by strictly positive functions P(a; s) for 
-1 < s < 1. 

On the other hand, positivity of the functions P(ct; s) 
for — 1 < s < 1 implies that the Glauber function, rep- 
resenting a formal limit P(a) = lim s _j.i P(a; s), is also 
non-negative and has singularities, not stronger than that 
of (5-function. Therefore, positivity of all the functions 
P(a; s) is a criterion for state classicality. 

With increase of the parameter s the function P(a; s) 
becomes more sensitive to state nonclassicality [for ex- 
ample, the Q function equal to P(a; —1) is always non- 
negative, but the Wigner function can take negative val- 
ues for certain states]. Therefore, the "order" of state 
nonclassicality (sensitivity of the observables to be used 
to detect the nonclassicality) can be characterized by the 
minimum values sq of the parameter s for which the 
phase-space function P(a; so) is not strictly positive (see 

Ref. my. 



s {p) = inf js I 3a : Tr (pS (a - a; s)) < oj . (16) 



For example, a single-photon state is extremely nonclas- 
sical: so(|l)(l|) = — 1 34^ - while a coherent state is a 
border between nonclassical states (with sq < 1) and 
classical states (formally with So > 1 — "nonclassicality" 
of classical state cannot be detected by any observable): 
s (\a){a\) = l. 

Fig. [5] shows the dependence of the nonclassicality or- 
der so of the stationary state of a single-atom laser on 
the system parameters. For v$ < the stationary state 
is nonclassical, with its nonclassical properties being de- 
termined mainly by the values of vq , similarly to predic- 
tions of the strong-coupling approximation (25j |. How- 
ever, in the region vq > the stationary state retains 
its nonclassicality, contrary to the characteristics of the 
approximate solution. The order sq of the nonclassicality 
almost does not depend on isq for uq > and is deter- 
mined mainly by cutting the " tails" of the photon number 
distribution. This type of nonclassicality corresponds to 
the inherent quantumness of single-atomic systems and 
arises for any parameters of the considered system. 

Fig. [S] illustrates the influence of the pump parame- 
ter eto on the stationary state nonclassicality. The char- 
acteristic properties of the nonclassical behavior resem- 
ble the predictions of approximate solutions (see, e.g., 
Refs. [23,133]) for different regimes of single-qubit laser 
operation. However, the stationary state remains non- 
classical even for "classical" regions. 
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FIG. 5: Stationary state nonclassicality order so(p/). Pa- 
rameters: Oo = 1, po = o-o + v o + 1; V = 5 (dashed line), 
50 (dot-dashed line). Solid line corresponds to the strong- 
coupling approximation. The values so = 1 corresponds to 
classical states, so < 1 — to nonclassical states. 
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FIG. 6: Stationary state nonclassicality order so(p/) (black 
lines) and correlation function f/ 2 ' (gray lines). Parameters: 
j/o = l (solid line), (dashed line), —0.5 (dot-dashed line); 
po = a,Q + vo, r) = 50. Regions of single-qubit laser operation: 
(a) — linear, (b) — quantum nonlinear, (c) — lasing, (d) — 
self-quenching, (e) — thermal. The inequality so < 1 (state 
nonclassicality criterion) holds for all the regimes, including 
the ones with > 1. 



VI. EFFECTIVE NONLINEAR TRANSITION 
PROBABILITIES 

As stated above, Eq. describes the balance be- 
tween energy dissipation from the system atom+field and 
pumping. Here we show that the second equation for de- 
termination of pn and P22 [Eq. ([3])] can be interpreted 
as the balance between the number of photons, absorbed 
from the field mode and emitted into it. 

To make the consideration more clear, we recall the 
scmiclassical description of an ordinary laser, consisting 
of a single mode and a large number of emitters. The 



average number of photons absorbed (N12) and emitted 
(N21) by each atom per unit time depends on the aver- 
aged number of photons (n) in the mode linearly and is 
determined by the constant Einstein coefficients 0: 



N12 = w(n)ipi, 



N21 = w{(n) 2 + l)p2, 



(17) 



(18) 



where p\ and P2 are the probabilities of finding the atom 
in the ground and excited states respectively; w is the 
spontaneous emission probability. The subscripts " 1" 
and "2" are used for (n) in order to take into account 
energy conservation, leading to change of the number 
of photons in the mode after absorption or emission: 
(n)i — (71)2 + 1- Then the steady state condition can 
be formulated as equality of the net number of photons 
emitted by atoms and the number of photons lost from 
the cavity: 

w((n) + l)(pa - pi) = 2/c«n) + l)(p x + p 2 ). (19) 

In the case of a single-qubit laser the quantities N12 
and N21 correspond to the transitions |n+l)|l) — > \n)\2) 
and vice versa, respectively. The probabilities of these 
states are equal to pn(n-\- 1) and p22{n). Therefore, the 
stationary state equation, analogous to Eq. (1191) . must 
have the following form [see Fig. H£b)]: 

(n + l)w n {p22(n) - pn{n + 1)} = 2k(h + l)p/(n + 1). 

(20) 

Comparing Eqs. © and (|2U1) . one can see that the tran- 
sition probability w n depends on the field intensity in the 
following way: 



K 



Mo + (n + 1) { 1 - 



(n + 2)p f {n + 2) 
(n + l)p f (n + l) 



For n <C Po the transition probability w n is approxi- 
mately constant, as it should be for ordinary spontaneous 
and induced transitions. However, for large photon num- 
bers n it becomes strongly intensity-dependent: 



Mo 



V 



V 



(21) 

and decreases with growth of n in such a way that the 
total transition probability (n+ l)w n tends to a constant 
value (Fig. [7]): (n + l)w n — > g 2 /n. Also the probability 
decreases with growth of the pumping rate TLyi — 4koq. 

Such single-atom blockade of intracavity photon emis- 
sion can be explained by the fixed " capacity" of the two- 
level atom, which is the only pumped object in the model 
considered. The atom can store only one excitation. This 
means that the system can accept only one energy quan- 
tum from the pump during the characteristic interaction 
time. Therefore, however large the probability of the in- 
duced transition of the atom from the excited state to 
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FIG. 7: Dependence of the effective total interaction-induced 
transition probability (n+ l)w n (expressed in terms oi g 2 /k) 
on the number of photons n, present in the mode: solid lines 
— numerical calculations; dashed line — approximate ana- 
lytical expression [see Eq. (|2ip ]. valid for large n; dot-dashed 
lines — total transition probability for w n = const. Parame- 
ters: black lines — vq = 1, a% = 1, fj,o = 3, r\ = 5 (coherent 
regime); gray lines — vq = 5, a 2 , — 5, fio = 200, n — 5 
(strong-dephasing regime) . 



the ground state can be, only one photon can be cre- 
ated in the mode during one such period. The system 
is effectively saturated by one photon, and the observed 
total transition probability is completely determined by 
the pumping, interaction and decay constants and does 
not depend on the field intensity. 

Mathematically, a coherent interaction between the 
atom and the mode leads to a correlated stationary state 
with the average photon number in the mode depend- 
ing on the state of the atom and, therefore, to effective 
suppression of interaction-induced transitions (the " net" 



transition probability decreases) . It should be noted that 
in the regime of strong atomic state dephasing (F" 3> K, 
Ho 3> 1) the correlation is rapidly broken, and the spon- 
taneous and induced transition probabilities behave in 
the ordinary way even for quite large n (see Fig. gray 
lines). 

VII. CONCLUSIONS 

To summarize, we have provided an analytical descrip- 
tion of the stationary state of a one-atom-one-mode sys- 
tem with incoherent pumping. The description captures 
both the features characteristic of the strong-coupling ap- 
proximate description and several new properties, such as 
stationary-state nonclassicality for all values of the inter- 
action parameters. The stationary state is shown to be a 
phase-averaged eigenstate of a special kind of deformed 
annihilation operator and, thus, to represent a phase- 
averaged nonlinear coherent state. The properties of the 
deformed annihilation operator and the obtained nonlin- 
ear coherent state are completely determined by the in- 
teraction parameters ao and Pq and the state-dependent 
operator d(n), diagonal in the Fock-state basis. The 
operator d(n) is constructed on the basis of an itera- 
tion scheme characterized by such important properties, 
as unconditional stability and independence of bound- 
ary conditions. Both numerical and uniformly applicable 
approximate analytical solutions are constructed on the 
basis of the iteration scheme. Interpretation of the sys- 
tem evolution equations in terms of spontaneous and in- 
duced transitions provided in our work reveals the inher- 
ent quantumness of a single-atom laser, which manifests 
itself in strong dependence of the transition probabilities 
on the field intensity and in a specific saturation effect. 
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